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We generalize the method of Van Hove so as to deal with the case of 
non-ordinary statistical mechanics, that being phenomena with no time-scale 
separation. We show that in the case of ordinary statistical mechanics, even if 
the adoption of the Van Hove method imposes randomness upon Hamiltonian 
dynamics, the resulting statistical process is described using normal calculus 
techniques. On the other hand, in the case where there is no time-scale 
■ separation, this generalized version of Van Hove's method not only imposes 

J> . randomness upon the microscopic dynamics, but it also transmits randomness 

y> : to the macroscopic level. As a result, the correct description of macroscopic 

\ dynamics has to be expressed in terms of the fractional calculus. 

PACS number(s): 05.40.+j, 05.45.+b, 05.60.+W 

OO : 
On . 

^ ! I. INTRODUCTION 

The physical paradigm of statistical physics is Brownian motion, which involves diffusion, 
dissipation and the fluctuation-dissipation relation tying the two together. The dynamical 
model of this process was provided by Langevin in 1908 using a stochastic differential equa- 
>- ! tion. In spite of this long history it seems apparent from the nature of randomness that such 

macroscopic stochastic equations are incompatible with the continuous and differentiable 
character of microscopic Hamiltonian dynamics. However, it is widely believed that Brow- 
nian motion can be rigorously derived from the totally deterministic Hamiltonian models of 
classical mechanics. Part of the reason for this conviction has to do with the wide use made 
in literature of Van Hove's method In one form or another, many of the attempts 

currently made to establish a unified view of mechanics and thermodynamics || can be 
traced back to the method of Van Hove. The result of this method depends on whether we 
adopt the Heisenberg perspective, corresponding to the time evolution of observables, or the 
Schrodinger perspective, corresponding to the time evolution of the Liouville density. In the 
former case, the usual outcome is the derivation from mechanics of an ordinary Langevin 
equation. In the latter case, the adoption of the Van Hove method yields master equations. 
We focus here especially on the conventional diffusion equation, with the diffusion process 
described by a second-order spatial derivative. 

In the Heisenberg perspective, after averaging over an ensemble of realizations of the 
stochastic force, the relaxation process is described by an exponential function. In the 
Schrodinger perspective, the mathematical representation of the diffusion process is given, as 
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we have said, by a second-order spatial derivative of a density function. Therefore, the math- 
ematical description rests on either ordinary analytical functions (exponential functions) 
describing the dynamics, or on conventional differential operators (second-order derivatives) 
describing the phase space evolution. This is probably the reason why there is no mention of 
the fact, and indeed no perception, that the Van Hove method is equivalent to introducing 
stochastic dynamics, namely discrete processes and unpredictable non-differentiable jumps, 
in the dominion of continuous and different iable Hamiltonian dynamics. 

Actually the differentiable nature of the macroscopic picture is, in a sense, a natural 
consequence of the microscopic randomness |5|,|B|], and of the related non-differentiability 
as well, due to the key role of the central limit theorem. Recall that in the central limit 
theorem the quantities being added together are statistically independent, or at most weakly 
dependent, in order for the theorem to be applicable and Gaussian statistics to emerge. Once 
a condition of time-scale separation is established, in the long-time limit the memory of the 
non-differentiable character of microscopic dynamics is lost, and Gaussian statistics result. 
This also means that use can be made again of ordinary differential calculations on the 
macroscopic scale, even if the microscopic dynamics are incompatible with the adoption of 
ordinary calculus methods. 

On the other hand, in the case where a time-scale separation between the macroscopic 
and the microscopic level of description does not exist, the non-differentiable nature of the 
microscopic dynamics is transmitted to the macroscopic level. An illuminating example is 
given by the paper of Ref. 0, which shows that a diffusion process generated by a fluctuation 
with no time-scale at the macroscopic level generates a diffusion process well described by 
a fractional Laplacian The paper of Ref. addresses also the intriguing problem of 

making compatible a Levy process, which has an infinite second moment, with the dynamical 
approach to diffusion. The dynamical approach to diffusion rests on steps of finite length, 
and consequently results in finite second moments. This problem has been addressed in a 
variety of ways JTU| , |TT|| , ranging from taking into account the finite size L of the sample, 



within which the Levy flight takes place [10], to the case where the probability density 



is truncated [[l ![. We rather follow the prescriptions of 0, which is a form of Levy walk, 
where the individual jumps are not instantaneous and involve a time cost, thereby making it 
possible for us to adopt a diffusion picture with finite moments. All this has been discussed 
in Ref. 0, and here the results of Ref. [|7] are made compatible with the infinite moments 
of the Levy statistics by means of the generalized Van Hove method. Another case where 
the adoption of a fractional calculus is made necessary concerns the generalization of the 
exponential form of relaxation according to the recent prescriptions of Nonnenmacher and 
co-workers |JT2"|,|r3|] , involving the concept of a fractional time derivative. The main aim of this 



paper is to prove that the fractional derivative, both the space fractional derivative of Ref. 
@ and the time fractional derivatives of Refs. Jl2],|l3j are made compatible with Hamiltonian 
deterministic dynamics by means of a generalized version of Van Hove's method. 

In fact, if we make the traditional classical assumption that the microscopic dynamics 
follow the Newton prescription, we also have to address the problem of how to make this 
prescription compatible with randomness. The solution to this problem requires that an ex- 
tension of the Van Hove method be discovered. This naturally yields the working hypothesis 
that the fractional derivatives, currently used to describe macroscopic transport processes 
P JT2]JTB| can be regarded as the macroscopic manifestation, in the absence of time-scale 
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separation, of either non-differentiable microscopic dynamics, an assumption that would 
violate the applicability of Hamiltonian dynamics to this domain, or of a Hamiltonian de- 
scription that loses differentiability through a kind of filtering described by a generalized 
Van Hove approximation. In other words, just as the Van Hove method makes Brownian 
motion compatible with Hamiltonian dynamics, a generalized Van Hove method is used 
here to make the macroscopic fractional calculus compatible with microscopic Hamiltonian 
dynamics. The main difference with the case of ordinary statistical mechanics is that in 
this non-ordinary case the non-differentiable nature of microscopic dynamics, either natural 
or forced by the adoption of the Van Hove prescription, is transmitted to the macroscopic 
level where it takes the shape of fractional derivatives. The main purpose of this paper is 
to substantiate with convincing arguments this working assumption. 

The outline of the paper is as follows. In Section II we express the Van Hove method in 
a form, equivalent to the ordinary method, but more convenient for the generalization that 
we plan to develop in this paper. In Section III we show that in the Heisenberg picture the 
generalization of the Van Hove method results in the fractional derivatives currently used 
by Nonnenmacher and co-workers to study polymer dynamics. Section IV and Appendix 1 
are devoted to the Schrodinger picture, and show that the Van Hove method in this picture 
leads to the fractional derivative introduced by Seshadri and West || (see also West and 
Grigolini ||) to describe Levy processes. Section V is devoted to concluding remarks. 



II. THE VAN HOVE LIMIT: AN EXAMPLE FROM A SPECIFIC MODEL 



As pointed out in the Introduction, the Van Hove limit [|TJ of a microscopic process 
turns out to be a key ingredient for the derivation of statistical mechanics from microscopic 
dynamics. We review this limit and how it works in the case of a specific model rather than 
in general. 

Let us start by recalling the meaning of this limit. Consider the integro-differential 
equation 

<m -A 2 [* K(t -t')f(t')dt', (1) 



dt Jo 

where / is some quantity of interest, if is a memory kernel and A is a parameter. Equation 
(H) is a typical non-Markovian equation obtained in studying physical systems coupled to 
an environment, and whose environmental degrees of freedom have been averaged over. In 
this case the parameter A can be regarded as the strength of the perturbation induced by 
the environment on the system of interest. 

In the literature, wide use is made of the Markov approximation which replaces the 
integro-differential equation (JJ) with the rate equation 



A / K(t')dt')f(t). (2) 



df(t) 
dt 

The Van Hove limit [Q] consists in making the limit A — > 0, t — > oo in such a way that the 
product X 2 t is kept constant. That is, setting x = X 2 t and F(x) = f(t), Eq.(|TJ) becomes: 

dF(x) /VA 2 



dx 



f K{t')F{x - \H')dt'. (3) 

J 
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Now, the adoption of the limit 



const = x = lim A t, (4) 
A^O 

t — ► oo 

makes it possible for us to replace the time convoluted form in Eq.(^) with 

df{t) 
dt 

where 



-A 2 r/(t), (5) 



roo 

T 



/ dtK(t). (6) 
Jo 



Eq.(^|) gives an exponential solution, the same as that obtained by means of the Markov 
approximation (0). 

Taking the limit A — > corresponds to assuming that the coupling of the system to the 
environment is weak, while the limit t — * oo means that the observation time is much larger 
than other temporal scales present in the system. Specifically this time must be larger than 
the microscopic time r. This remark allows us to reformulate the Van Hove limit in a slightly 
different way, more suitable for our purposes. First, instead of taking the t — > oo limit, we 
shall take the limit r — * 0. Also, we shall replace the limit A — > with the equivalent limit 
V — > oo, where V is a coupling constant, to be specified in the following. The quantity to 
be kept constant in carrying out the limit is just the product V 2 t. Notice that connecting 
V, t and A as A = Vr makes it possible to keep V <C r _1 (in such a way that A < 1) and 
at the same time to set V 2 — > oo so as to ensure V 2 t — > const. 

The meaning of Eqs. (|3|) to (|5|) can be illustrated by adopting this equivalent perspective 
and can be shown to be a way of disregarding the time evolution of the system at both very 
short and very long times, where the deviations from the exponential relaxation show up 



1 15]. In the Hamiltonian case the relaxation cannot be rigorously exponential [0 thereby 
preventing a satisfactory connection between microscopic dynamics and stochastic physics 
[T?| . We now show this benefic effect of the Van Hove method with the help of an illustrative 



example, where the origin of Brownian motion stems from the average over an initial statis- 
tical distibution rather than from chaos 0. Consider a chain of 2N + 1 linear harmonic 
oscillators all with equal spring constants and described by the Hamiltonian 

N 2 u N 

"=.E£ + 5.£>'-*>' (7) 

i=—N 1 i=—N 



where m 8 = m if i ^ and mo = M. Vitali and Grigolini [IS] prove that the correlation 
function for the momentum p of the system of interest 

*.M - ^ (8) 
satisfies the following integral equation 
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where 



$o = -A 2 / ^(t-O^o^X 
Jo 



A 2 = 2k/M 



(9) 



(10) 



and $i(i) represents the correlation function of the stochastic force in the corresponding 
generalized Langevin equation (see Section III). The parameter Ai plays the same role as 
that of the earlier mentioned coupling constant V, while r is given in this specific model by 
the expression: 



In the long-time region, the solution to Eq.(|9[) is given by |18| 



(11) 



*o(t) 



1-2/i 

\ — 1± 
l-2fj, 



exp 



exp 




2\i r°° sm(xu! t)\^' 

7T Jo ( 



2/x)x 2 + fl 2 



-dx 



(i-#n^) sin (" 0( "i 



where 



and 



/i = m/M 



4k/m. 



(12) 



(13) 



(14) 



The asymptotic approximation (|T2J) shows that the momentum autocorrelation function 
consists of the sum of an exponentially decaying term and a non-exponential oscillatory 
contribution that decays as an inverse power law in time. Notice that ujq is a Debye-like 
cutoff frequency that is related to the microscopic time-scale r as r = 2/uq. The general 
form of this solution had also been observed and discussed by Zwanzig || 

Now, the Van Hove limit expressed in terms of the parameters of the model becomes: 
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A? 



lim 

T -> 

2 



A^r = lim cuo/x. 
fi — ► 

OO UJq — > oo 



(15) 



Notice that this limit can be realized assuming m —> and M — > as m 1 / 2 . The closer 
to zero the mass m, the better the physical condition M > m is fulfilled (the better the 
macroscopic description of the system of interest). Consequently, the oscillatory tails in 
Eq.(|l2"D cancel and an exact exponential relaxation is recovered: 



$ (t) = exp(-7t). 



(16) 



The rationale for this result stays in the formal similarity of Eq.(^) with Eq.(|]). Letting, 
as mentioned before, A = Air, it is a simple matter to verify that A — > as m 1 / 4 and 
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the ordinary Van Hove limit (Q) can be applied. More in general, we can assume that the 
relation fllEf) holds true independently of the model, and we can regard the parameters \x 
and Uq as free parameters of the theory. If we adopt this view, the exponential decay is 
recovered by making the Van Hove limit as described by the r.h.s. of Eq . flT5|) . 

In the following, we shall refer to the Van Hove limit, rather than using the known version 
of Eq.(||), as expressed in our non-conventional form: 

7 = lim A 2 r = Limv H A 2 T. (17) 
r -> 

A 2 ^ oo 

A final remark concerns the fact that the same result, Eq.(|l^), could have been obtained 
by applying the Markov approximation to Eq.@. In that case, the basic assumption would 
have been an infinite time-scale separation between the microscopic time-scale r and the 
macroscopic scale defined as the inverse of the frequency Q 2 = 2k/M, that is: 




Notice that the time-scale separation between the system of interest and bath is rendered 
infinitely large exactly by the same limits as those used in order to carry out the Van Hove 
limit: m — > and M — > as m 1 / 2 . This demonstrates that the Van Hove limit and the 
Markov approximation are essentially equivalent to one another. 



III. THE HEISENBERG PICTURE: ORDINARY AND FRACTIONAL 

RELAXATION 



It is well known |19| that the generalized Langevin equation 

v = -A 2 [ - t')v{t')dt' + f{t) (19) 
J o 

corresponds to the following hierarchy of correlation functions: 

= -A 2 +1 f $ i+1 {t - t')^(t')dt'. (20) 
Jo 

Eq.(|) is the i = case of Eq.([2g). 

In order to derive both normal and anomalous relaxation properties, we are interested 
in making a non-trivial choice of the correlation function $i(£). However we also want 
our choice to be compatible with a completely dynamical approach. Therefore we need to 
identify the conditions necessary to assure that both these constraints are satisfied. First, 
by means of a Laplace Transforms and the convolution form of Eq. (|20|) , it is easy to prove 
that $i can be represented in the form of a continued fraction: 

= ^2 > (21) 

z + — 2 



z + 



z + ... 



where $1(2) denotes the Laplace Transform of This structure is valid for $i(0) = 1 

and this is a first requirement to fulfill. Then, we recall that the expansion parameters Af 
can be expressed in terms of the moments 

Sn- {m) (0), (22) 

where L is the Liouvillian operator driving the time evolution of the Liouville density and 



is the first state in the Mori chain [|T{J. An elegant way of expressing the parameters A 
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in terms of the moments s n has been established by Grigolini et al. p0[ . This implies that 
$i(i) must be infinitely different iable. Finally, the symmetry properties of the Liouvillian 
imply that also the condition 

S2n-l = (23) 

has to be fulfilled. 

Therefore, we decide to focus our attention on the choice: 

^ = (T* + (24) 
Note that the moments s n are nothing but the n-th order time derivatives of $i(t). Thus, 



it is straightforward to prove via successive time differentiation of (24]) that the odd mo- 
ments vanish. This is an important mathematical property necessary to make the relaxation 
process compatible with Hamiltonian dynamics. An important example of relaxation incom- 
patible with Hamiltonian dynamics is in fact the exponential relaxation: A case exhaustively 



discusses by Lee in his brilliant 1983 paper [21 1. This paper proves therefore the importance 



of fulfilling the constraints of Eq. fl23|) for a generalized Langevin equation to be compatible 
with Hamiltonian dynamics. This condition is here fulfilled by adopting the choice of Eq. 



A. Ordinary statistical mechanics 

Ordinary statistical mechanics can be recovered from the generalized Langevin equation 
(fOp, using Eq.(|24]) along with the integrability condition on the power-law index, 

P > 1. (25) 

The microscopic time-scale r is given in this case in terms of the parameter T, 



T -, £ T ' ^ r (^) T (26) 

- ' "\r' + p)m - 2 r(f) 

Therefore, the Van Hove limit in the form ([H]) is achieved in the limit T — ► 0, 

7 = ^£l^I L im v „A;r. (27) 



r(f) 
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This limiting procedure results in exponential relaxation for the correlation function $q(£) 
and allows us to safely interpret Eq.(19) as identical to the ordinary Langevin equation, 



v = ->yv + f(t). 

This is the traditional result obtained using the Van Hove method. 



(28) 



B. Non-ordinary statistical mechanics 

In the case where the power-law index of the correlation function is in the interval 

< j3 < 1, (29) 



the non-integrability of the correlation function ( p^ ) prevents us from adopting the above 
approach. In the case < (5 < 1, we are forced to look for a different procedure to go from 
a microscopic to a macroscpic description of the system. This procedure can be derived in 
a natural way from the original Van Hove limit. Let us consider the following limit, 

Q = lim A\T? = Lim GVH A?T^. (30) 
T -> 

A\^ oo 

We shall refer to this limit as to the generalized Van Hove limit. 

Adopting the ansatz (|30|), and inserting Eq.(^4j) into Eq. (p~9|) , we obtain 

*„ = -Lim GV „A;ro£ (T2 + {t \ t , )2)m M^df 

-- Q [jTwW- < 31) 

For dimensional reasons, it is convenient to write 

Q = V 2 t?. (32) 

Notice that in general the correlation function is related to the waiting time distribution 
ip{t) of the process under study as 

$ (t) = 1 - f M)dt' (33) 
Jo 

and therefore Eq.([H]) can be rewritten: 

^ = Q ll jrhy Mt ' )dt '- (34) 

Now, we want to compare this last expression with the fractional relaxation equation ob- 
tained by Glokle and Nonnenmacher [12]. The fractional relaxation equation for a function 
$ml(£; v) is given by 
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1 dr v 

$ml(£; v) - $ ML (0; v) = — — -$ml(*; v) (35) 

at v 



where the symbol denotes the fractional integral (see Appendix 1), 



dt-^ J v y r(z/) 7 (t -t') 1 -"' 



The solution of Eq.(Rq) is known [13] and is given by the so called Mittag-Leffler function, 



"> = | rpTi) (;)"*■ (37) 

which exhibits stretched exponential behavior at short times and inverse power-law relax- 
ation at long times. 

Setting (3 = 1 — v, it becomes possible to compare the fractional relaxation equation with 
Eq.(||) to obtain 

m = -QV(1 - (3)r l -P [$ ML (t; 1-/3)- $ ML (0; 1 - (3)} . (38) 

where the superscript 1 — (3 corresponds to the value for v to insert into the form of the 
solution QUI B|) . Recalling also Eq . fl32|) , we finally obtain for the waiting time distribution 
function 

m = -(V 2 r)T(l - P) [$ ML (t; 1-0)- $ ML (0; 1 - (3)} , (39) 

which mantains the interesting properties of the Mittag-Leffler function earlier pointed out. 

It is important to stress that the resulting analytical function has been used by Glockle 
and Nonnenmacher |22[] to fit with a very remarkable accuracy the relaxation curves of stress 
experiments on glassy material. This suggests that the dynamical randomness without time- 
scale separation takes the shape of a time fractional derivative and becomes experimentally 
detectable at the macroscopic level. 

IV. SCHRODINGER PICTURE: GAUSSIAN AND LEVY DIFFUSION 

Let us now consider the equation of motion for the one- dimensional probability density 
p(x,t): 

§- t P(x,t) = {e)e q f Q ^{t')^p{x,t-t')dt'. (40) 
This equation refers to the process: 

x = £. (41) 

and $^(t) is the autocorrelation function of the £— fluctuations. The integro-differential 
equation of motion for the probability density is obtained from the Langevin equation by 
adopting the Zwanzig projection operator method in the form discussed by Grigolini |23| 
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and it is exact under the following two conditions: (i) the dynamics of £ is independent of 
that of x and (ii) the system producing the £— fluctuations is a two-state system. 

Notice that (ii) does not necessarily mean that the variable £ is dichotomous. The case of 



anomalous diffusion generated by intermittent maps |24| is an illuminating example, where 
the change from a continuous two-state fluctuating variable to the dichotomous case does 
not produce significant effects on diffusion. From the point of view of our generalized Van 
Hove limit, however, the difference in the two conditions is essential. As we shall see in the 
next subsection, the replacement of £ with a dichotomous variable with the values W and 
—W is an essential ingredient of the generalized Van Hove method. 

For reasons which will shortly become clear, it is convenient to write the probability 
distribution at time t — t' in terms of the probability distribution at time t as 

/oo 
F(x - x', -t')p(x', t)dx', (42) 
-oo 

where the function F is a propagator. Equation([R]), with the non-trivial choice (|24]), be- 
comes 

d r°° ft T 13 d 2 

-p(x,t) = (ZXilJx'l dt' {T2 + t , 2)p/2 ^F(x - x>, -t')p{x\ t). (43) 

We shall discuss both the case of ordinary diffusion, f3 > 1, and the case of anomalous 
diffusion, < j3 < 1. 

A. Ordinary statistical mechanics 

Let us assume that the power-law index in the autocorrelation function is 

P > 1 (44) 
and let us make the choice of a time-independent propagator 

F(x-x',-t') = 5(x-x'). (45) 



Equation(f43D then becomes 



m p(z, t) = (eU {I *t {T2 + f y /2 ) gj**, 0. (46) 

The microscopic time-scale, r, of the autocorrelation function, $^(t), present in Eq.(^) is 
finite for f3 > 1 and we can apply again the conventional Van Hove limit in the form 

D = lim (£ 2 ) e(? T. (47) 
T -> 

d 2 ) eq -> oo 

Substituting this expression into Eq.(^) we obtain the standard diffusion equation 

§- t Pix ' t)=D & P ^ t) - (48) 
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Note that the case where the correlation function is exponential, although in conflict 
with the Hamiltonian constraints fT6|, assists us in further interpreting the meaning of 
the Van Hove limit. As shown in Appendix 2, using an exponential correlation function 
in Eq.(f43D along with Eq. (^5|) , reduces Eq . ([43D to the telegraphic equation but ordinary 
diffusion can be still recovered by making use of the Van Hove limit. 



B. Non-ordinary statistical mechanics 

Let us consider now the case where the power-law index is in the interval 

< p < 1. (49) 

In this case, we define the generalized version of the Van Hove procedure as the recipe leading 
to the largest component of the diffusion process. Therefore, as a first step, we make the 
variable £ dichotomous, by assigning to it the values W and — W with the basic condition 
that the actual values of the true variable are included in the interval [-W, W] (see Fig. 1). 




-w w \ -w w \ 

FIG. 1. Sketch of how the distribution P(£) of the £ variable is modified so as to make the 

variable £ fully equivalent to a dichotomous variable. 

This leads us to assume for the form of the propagator F, 

F(x - x', -t') = 5{Wt' -\x- x'\), (50) 
which is clearly time-dependent. Equation ( |43|) then becomes 

d . . W p+l TP r°° , ,,, ^[W 2 T 2 +\x-x'\ 2 ]-(3((3 + 2)\x-x'\ 2 

It is evident that Eq.(|5l|) differs from Eq.(f|3|) by a correction term that disappears in the 
generalized Van Hove limit. The generalized Van Hove limit is defined in this case as: 

Q = lim T P W 1+P (52) 
T -> 

W -> oo 

and noticing that this implies W 2 T 2 — > as W~ 2 ^, Eq.(|5lD can be rewritten 

m pM = —^ Q L dx ix-^xW 2 - (53) 
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This form of diffusion equation coincides with the West-Seshadri equation |7||| for a centro- 
symmetric Levy process. The fact that Eq.(|53[) allows for a solution in the form of a Levy 
process becomes apparent by taking its Fourier transform. Letting 

11 = 0+1 => 1< fj, < 2 (54) 

and using |]25|] , we obtain for the Fourier transform of an inverse power law 

2 r(-/i)|fcrcosf-^V (55) 




7T V r/ ' ' V 2 

Therefore, making use of the Convolution Theorem for Fourier Transforms, Eq.fl5"3|) becomes: 

^<P(k,t) = -b\k\»<P(k,t). (56) 

Here, 4>(k,t) is the Fourier Transform of p(x,t), that is, the characteristic function, and the 
parameter b is given by 



V Z7T 



COS 



~2~ 



> 0. (57) 



The solution to Eq. (|56|) , with the initial condition (p(k,t = 0) = 1, necessary for the inverse 
Fourier transform to be defined as a probability density, is written as 

(f,(k,t) = e - 6|fc|Mi , (58) 

which indeed corresponds to the definition of the characteristic function for a Levy process. 
Equation (|5^) can also be cast in the form of a fractional differential equation U 

^p{x,t) = {-i)»bVMp(x,t), (59) 

where T>^> is a proper definition of fractional derivative (see Appendix 1). We see, in con- 
clusion, that with the adoption of the generalized Van Hove method the standard diffusion 
equation yielded by the ordinary Van Hove method is expressed in terms of a fractional 
derivative. More in general, we expect that the ordinary Fokker-Planck equation can be 
replaced by a generalized expression resting on fractional derivatives ||. 

In the next Section we shall give more support to our conviction that the generalized 
form of Van Hove method, here adopted to derive the spatial version of fractional derivatives 
of (|59|) is intimately related to the Levy-Gnedenko theorem . We limit ourselves to point 
out that the adoption of fractional calculus to deal with processes of anomalous diffusion 



is becoming more and more popular, and we quote for the interested reader Refs. p7| , |28 
whose results, however, must be compared to the conclusion of Ref. as well as to those 
of this Section with some caution |50] 
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V. CONCLUDING REMARKS 



The adoption of the Van Hove limit is essentially a sophisticated way of making the 
Markov approximation. The Markov approximation, in turn, establishes the physical condi- 
tion necessary to make Hamiltonian dynamics compatible with stochastic physics. However, 
when this method is applied to microscopic dynamics to derive ordinary statistical me- 
chanics, there is no clear perception of establishing dynamic properties inconsistent with 
Hamiltonian dynamics. This is so because even if the correlation functions are made expo- 
nential by forcing the Markov approximation into microscopic dynamics, so as so to become 
incompatible with both classical [JIB and quantum |0J mechanics, they are still differentiable 
functions. 

From the point of view of a single trajectory, the realization of the Brownian condition 
implies, from a rigorous mathematical point of view, the breakdown of the condition of differ- 
entiability. However, even in this case the perception of a conflict with Hamiltonian dynamics 
is blurred by the adoption of a statistical perspective. Within the Schrodinger-like picture, 
namely the picture where we observe the motion of an ensemble of trajectories rather than 
that of a single trajectory, ordinary diffusion is produced by a dynamical operator propor- 
tional to a second-order spatial derivative. Again, a condition of ordinary differentiability is 
ensured! In conclusion, the existence of a time-scale separation between microscopic dynam- 
ics and the experimental observation, usually made at the macroscopic level, ensures that 
the results of the observation process can be predicted by means of theoretical prescriptions 
based on ordinary mathematical procedures resting on the differentiability assumption. 

A totally different condition is generated if the time-scale separation is not adequate. 
We have separately discussed the Heisenberg and the Schrodinger picture. From the results 
of Section III, devoted to the Heisenberg picture, we see that when the condition of time- 
scale separation is insured, the ordinary Van Hove method can be applied, and standard 
exponential relaxation follows. In the absence of the condition of time-scale separation the 
method of Van Hove must be generalized and this, in turn, yields a generalization of the 
exponential relaxation, a condition that according to Glockle and Nonnenmacher |12,13| 



turns out to be very efficient for describing such non-standard physical processes as polymer 
dynamics. 

From within the Schrodinger picture, discussed in Section IV, the lack of a time-scale 
separation and the consequent requirement for a generalization of the Van Hove method, 
yields the striking replacement of the ordinary Laplacian operator with the fractional Lapla- 
cian discussed in Appendix 1. We have to point out that even in this case randomness has 
a very subtle origin, implying a departure from the differentiable condition of Hamiltonian 
dynamics, and the replacement of the continuous variable responsible for microscopic fluctu- 
ations with a dichotomous variable. This process, changing a two-state physical dynamics, 
compatible with a Hamiltonian picture, into a non-differentiable process is clearly described 
by the sketch of Fig. 1. We firmly restate here that the replacement of a continuous variable 
picture with a discrete representation is done also in the case of ordinary statistical mechan- 
ics, where this way of forcing randomness within differentiable dynamics does not imply a 
departure from differentiable dynamics at a statistical level. When the condition of time- 
scale separation is not available, the Van Hove procedure of forcing microscopic dynamics 
to become random has a striking manifestation at the macroscopic level under the form of 
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the fractional Laplacian of Appendix 1. 

To stress the main conclusion of the paper, we set aside the non-Hamiltonian cases, and 
especially those where stochastic processes are assumed to be already at work, whatever their 
origin might be. This is where the microscopic dynamics are genuinely not differentiable, 
and whether or not this lack of microscopic differentiability has macroscopic manifestations 
depends on whether or not this non-differentiability can be transmitted to the macroscopic 
level as a result of memory The present paper sheds light also on this interesting, but 
less fundamental issue. A much more important condition is that where the dynamics are 
described by a Hamiltonian. In this case the microscopic dynamics are differentiable, and in 
principle, there would be no microscopic randomness to transmit to the macroscopic level. 

However, one of the main tenets of the current literature on the derivation of statistical 
mechanics from dynamics rests on the conviction that randomness can be generated from 
hamiltonian dynamics either as a consequence of the action of infinitely many degrees of 
freedom or as a consequence of chaos. The Hamiltonian of Eq.([7|) is an example of the first 
kind, and the sense of the Van Hove limit in that case is that the very short and the very 
large time-scales are beyond the range of observation. Thus, randomness is a consequence 
of the observer's limitations. The physical condition illustrated by Fig. 1, on the contrary, 
might be related in some way to the case of chaos, in a sense that it is convenient to properly 
discuss here. 

Let us consider the well known case of a map resulting in intermittency for instance, 
that introduced by Geisel and Thomae ||32|| . The motion under study is characterized by a 
long permanence in conditions of regular motion with bursts of chaotic dynamics concerning 
the transition from one to the other laminar region. This is not a genuinely Hamiltonian 
dynamics and the generator of intermittency is in fact a map. According to the general 
discussion of PTJ, similar properties can also be exhibited by genuinely Hamiltonian systems 



such as the attractive model of the egg-carton two-dimensional potential investigated by 
Geisel, Zacherl and Radons [Q. Another interesting model of the same type is the three- 
dimensional Hamiltonian flow more recently studied by Zaslavsky, Steven and Weitzner 
34]. This paper, as well as that of Ref. p3| , can be regarded as an example of Hamiltonian 



derivation of the Levy processes and consequently, according to the point of view adopted 
in this paper, an example of a dynamical process described at the macroscopic level by a 
fractional derivative in spite of its Hamiltonian, and consequently, differentiable nature. 

It has to be stressed, however, that the theory behind this macroscopic derivation should 
be applied to the numerical treatment of these processes: This theory does not have any- 
thing to do with the ideally exact solution of a dynamical model resting on a continuous 
treatment and the assumption of differentiability at any order. The numerical treatment is 
characterized by round-off errors and, more importantly, by a discrete time representation, 
which forces the system to depart from the conditions ideally established by its Hamiltonian 
property. The intermittency of these systems are sufficient to create at least temporary 
conditions of non-differentiability which are then transmitted to the macroscopic level and 
changed, in the way described in this paper, into fractional derivatives. The generalized 
version of the Van Hove method serves the basic purpose of introducing conditions of micro- 
scopic non-differentiability without leaving the theoretical treatment and without entering 
the level of the numerical solution. In the case of the dynamical approach to Levy processes, 
the meaning of the generalized Van Hove method is closely related to the Levy-Gnedenko 
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generalized central limit theorem |pq| . This is so because forcing the fluctuating variable 
to become rigorously dichotomous imposes on the system enough coarse graining as to pro- 
duce randomness. We note also that in the case 1 < (3 < 2, T is proportional to the mean 
time duration the map spends in one laminar region. Consequently if t 3> T, the number 
N = t/T, which is very large, corresponds to the number of uncorrelated space transitions 
of intensity \x — x'\, with probability l/\x — x'\ 1+l3 , made by the system. The Levy-Gnedenko 
theorem insures that for N tending to infinity a Levy diffusion process is generated. From 
a physical point of view this has the same effect as forcing T to zero. 

We think that at this stage we are also in the right position to establish an appealing 
connection with exciting results of the research work of Refs. |5|,§]. Although some doubts 
are expressed by the authors of Ref. || on the role of chaos to generate Brownian motion 
(a role which might be sufficient but not necessary, as, in a sense, it is also shown by our 
chain model of Section II) we are inclined to believe that the connection established by the 
authors of Ref. || between Brownian motion and chaos is very attractive. However, it applies 
to the condition of time-scale separation which is expected to generate ordinary statistical 
mechanics. This paper is devoted, on the contrary, to studying dynamical cases where this 
time-scale separation is missing, and consequently, non ordinary statistical mechanics is 
generated. This might generate the wrong impression that the method of analysis adopted 
in Ref. |J, based on the use of the Kolmogorov- Sinai entropy [[35], cannot be applyed to 



the dynamical systems with no time-scale separation. Actually, it has been recently shown 
| 3~6fl that the Kolmogorov- Sinai entropy can be generalized so as to be made efficient also 



in the case of fractal dynamics. Furthermore, it was also recently shown |37| that the same 
entropic arguments naturally lead to the same form of spatial fractional derivative as that 
of Eq. (|59|) . This means that the Van Hove generalized method discussed in this paper is 



expected to establish a natural bridge between the non-extensive entropy of Tsallis and 
the fractional derivative of Eq. (|59|) much in the same way as the ordinary Van Hove method 
makes Hamiltonian dynamics compatible with the standard diffusion equation of Eq. (fiSD . 
and so with ordinary extensive thermodynamics behind it [ |39|| . 
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APPENDIX 1 

We want to define the integrals and derivatives used in the fractional calculus introduced 
in the text. First of all, let us recall the Riemann-Liouville definition of fractional integral. 
Let us assume (3 > and, following [lOfl, let us define the /^-fractional integral of the function 
f(x) to be 



dx~P Jy ' T{(3)Jc {x-yf-P 
We also define the j3 fractional derivative of the function f(x) as 
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dx$ Yin- (5) dx 11 Jc (x - j/)/*-^ 1 

where n is the smallest integer larger than (3, that is n = \0\ + 1. The constant c in the limits 
of the integrals is usually set to (Riemann definition) or to — oo (Liouville definition). It 
is easy to show that that for (3 integer both definitions reduce to the ordinary definitions of 
derivative and integral. 

An equivalent definition makes use of Fourier Transform |^TJ. Let us consider a function 
f(x) with Fourier Transform f(k): 

f(k) = -= / f(x)e ik *dx, (62) 



2,7V 



oo 



1 f°° f 

f(x) = -= / f(x)e- ikx dk. (63) 

V J —oo 

The ra-th derivative of f(x) can be written as 

D^f{x)=F- l [{-ik) n f{h)-x\ (64) 
and a possible way of generalizing this expression to the (3 derivative of f(x) is 

D^f{x)=T- l [{-ikYf{k) ] x\. (65) 
Eq.(|65D is equivalent to the convolution product: 

D^f(x) = -L r dW)(a; - y)/(y)dy, (66) 

VZ7T J —oo 



/2tt 

where 



d(fl(z) = -= / {-ikfe~ lkx dk. (67) 
v27r J —oo 



It is possible to calculate an explicit representation of the integral (p7|) and the result is []4T 



r(n - /3) cfe n J-oo (x - yf 



D W f( x ) = L —*L. r fiy)dy (69) 

J[ } r(n - /?) dx" L {x - y)/ 3 ^ 1 ' 1 ' 

where again n = [f3] + 1 and the +(— ) corresponds to evaluating the integral fl6"7|) in the 
upper (lower) complex k plane. It is apparent that Eq . (|68|) coincides with the Liouville 
definition of fractional derivative: 

Dff{x) = ^f(x). (70) 

A different way of generalizing Eq.flBip is the following. Let us start from the definition of 
second derivative 
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D^f(x)=^- 1 [(-k 2 )f(k);x} (71) 

and let us generalize it as: 

V^f(x)=^- 1 [(-k 2 ff(k)-x]. (72) 
Eq.(|72"D is equivalent to the convolution product: 

V^f(x) = ^=T c<M{x - y)f(y)dy (73) 

where 

c W{x) = L^L / \k\ 2f3 e- ikx dk. (74) 

If 

< /3 < 1, (75) 
the integral (|74]) can be evaluated explicitly |25 



C W( X ) = ( 1L (76) 

y^r(-2/3)cos(/37r)M 2/3+1 ' 



and setting 2(3 = p, (=>• < /U < 2), we get: 

_1W2 

-c M (a;). (77) 



1 2& (-ir /2 



|a^ +1 fi(fj,-l)Q 

Using now Eq. fl54|) , Eq.(|73|) and Eq. (|77|) , Eq.(|53"D can be expressed in terms of the fractional 
derivative T>^\ This gives Eq . (|59"D . Notice that Eq.(|59|) is consistent with Eq. fl56|) since 

T[c^\x)] k] = {-k 2 Y' 2 = i^kl". (78) 



APPENDIX 2 

Our aim is now to show that the conventional Van Hove limit produces in the case 
of exponential relaxation for the correlation function the same result as the traditional 
approach. Let us consider the following equation of motion for the probability density 
p(x,t): 

^P(x,t) = (v 2 ) eq J*$ v (t-t')-^p(x,t')dt f . (79) 

Eq. ([79"D is the same as Eq. (|40|) and the bath variable, now denoted with v, is not dichotomous. 
Assuming 
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-7* 



1 

7 = - 

T 



and differentiating both sides of Eq.(|79|), we obtain 



d 2 , 9 , „ . , 2v d 2 

eq dx 2 



—p(x,t) = -7— p(x,t) + (w 2 ) 



p(x,t), 



whose Fourier Transform is given by 



J^p(fc, *) + 7 J^(fc, t) + (v 2 ) eq k 2 p{k, t) = 0. 



(80) 



(81) 



(82) 



Both Eq.flSTD and (jS!^) correspond to the telegraphic equation. Considering the condition 



j)(M) = o, 



the exact solution of (|82|) is written as 

p(fc,t) = A 



Q! + (fc)i a_ (fc)t 

a_(A;) 



with A to be specified according to the normalization condition and a± given by 

n- -- ± 



2 2 



7 



Assuming now 



4(^%fc 2 
7 2 



< 1 



(83) 



(84) 



(85) 



(86) 



the solution 

p(k,t)=A 



can be expanded in a Taylor series, giving 
(v 2 ) eg k 2 \ 



cxp 



7 



7 - (v 2 ) eq k 2 /^ 



exp 



7 



7 



(87) 



Note that setting condition (|S6|) is compatible with the Van Hove limit resting on setting 
7 —>■ oo. At this stage we have to apply to Eq. ( j87l) the Van Hove limit, which in addition 
to setting 7 — > 00 rests also on making the limit (v 2 ) eq — > 00 in such a way that 



D = Lim 



VH" 



<f 2 > 



eq 



7 



We note that this procedure makes the terms O 
ordinary diffusion equation solution: 

p(k,t) = Ae~ mH 



(v 2 ) 2 eg 



(88) 



disappear so as to recover the 



(89) 



On the other hand, the same result can be obtained by applying the Van Hove limit directly 
to the equation of motion ( |79| ) with the assumption (|80| ) and also making use of both Eq.(f4"2"|) 
and Eq.([45|). This proves that ordinary Brownian diffusion rests on the Van Hove limit. 
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